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the skutterudite compound CeFe4Pi2 by Electron Spin Resonance. 
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Despite extensive research on the skutterudites for the last decade, their electric crystalline field 
ground state is still a matter of controversy. We show that Electron Spin Resonance (ESR) measure- 
ments can determine the full set of crystal field parameters (CFPs) for the cubic symmetry (/m3) 
of the Cei-2;Ra:Fe4Pi2 (R = Dy, Er, Yb, x < 0.003) skutterudite compounds. From the analysis 
of the ESR data the three CFPs, B^, Bg and Bg were determined for each of these rare-earths at 
the Ce''"*' site. The field and temperature dependence of the measured magnetization for the doped 
crystals are in excellent agreement with the one predicted by the CFPs derived from ESR. 



INTRODUCTION 

The filled skutterudite compounds RT4X12, where R 
is a rare-earth or actinide, T is a transition metal (Fe, 
Ru, Os) and X is a pnictogen (P, As, Sb) crystallize in 
the LaFe4Pi2 structure with space group /m3 and local 
point symmetry T^ for the R ions. Also recently a new 
skutterudite family, (Sr,Ba)Pt4Gei2, was found. [l| The R 
ion is surrounded by eight transition metal ions forming 
a cube, and twelve pnictogen ions that form a slightly 
deformed icosahedron. 2] These materials exhibit a broad 
range of strongly correlated electron phenomena. [3, 0, 
0] In addition, the antimonite members, are potential 
thermoelectric materials due to their enhanced Seebeck 
coefficient. [1,01 

It has been assumed for a long time that the descrip- 
tion of the electric crystalline field (CF) of the cubic point 
groups, T, T/i, O, Td, and Oh is the same for all of them. 
Recently, Takegahara et al., 8] studied the CF for cu- 
bic point groups using group theory and a simple point 
charge model and found that the above was not correct. 
Takegahara et al. noticed that due to the absence of two 
symmetry operations in the T and T/j groups, namely the 
C4 and C'2 rotations, the CF Hamiltonian (Hcf) al- 
lows for additional sixth order terms with an extra crystal 
field parameter (CFP), _Bg. Therefore, for Th symmetry, 
in terms of the Steven's operators[lOj Hcf should be 
written as 
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where the last term is absent in the ordinary cubic sym- 
metry Oh- Its presence does not affect the degeneracy of 
each sublevel when compared with that of the Oh group, 
but some eigenfunctions and eigenvalues may be appre- 
ciably different, jsj The knowledge of the CF levels, es- 
pecially the ground state, is essential to understand the 



role of the 4/-electrons in these compounds. However, in 
spite of the large amount of work invested, the C F g round 
state is still unclear in several of these systems. [Ill Il2| 

Electron spin resonance (ESR) has been used for 
more than half a century to examine a wide variety of 
compounds. 10] It is a very useful and highly sensitive 
technique to study spin correlations. It provides infor- 
mation about CF effects, site symmetry, valence of the 
paramagnetic ions, g-value, fine and hyperfine parame- 
ters, etc. Besides, the sample size required for ESR is 
typically less than ^ 4 mm^, i.e., much smaller than that 
needed for most other techniques. When the compound is 
not paramagnetic, ESR can still provide useful informa- 
tion by doping the matrix with a small amount of para- 
magnetic ions such as Nd^+, Gd^+, Dy^+, Er^+, Yb^+, 
etc. The ESR spectra of the impurities allow not only to 
learn about the impurity, but also to study the proper- 
ties of the host lattice. In cases where the first excited 
state is separated from the ground state by an energy of 
the order of the temperature at which the data is taken, 
a field induced change of the g- value [is'l and an expo- 
nentially activated T-dependence of the linewidth 14l.ll5| 
may be expected. Moreover, the ESR of an excited state 
could also be observed. ^16*1 Thus, by measuring the ESR 
at different frequencies and temperatures of various R 
impurities, one may obtain an accurate determination 
of their ground state and, in some cases, the full set of 
CFPs determining the overall splitting of the ground J- 
multiplet. 

Mesquita et al.., 17 1 and Martins et a/.,[l3| measured 
the ESR spectra of Cei_::,R^Fe4Pi2 (R = Nd, Dy, Er, 
Yb; X < 0.005) up to 4.2 K. Our data, taken in the same 
range of T, agree with those published previously. The 
data in Refs. [l7l.[l8j was analyzed assuming Hcf for the 
cubic group, i.e., Eq. |T]) without the i?g-term. In partic- 
ular, the unexpected g-vahxe of 6.408(3) measured for the 
Kramers doublet ground state of Er'^+ in CeFe4Pi2 can- 
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not be explained if the term 3^(0^ — Og) ^'^^ included 
in HcF- By using the Hcf given in Eq. ([T]) and measur- 
ing up to T = 50 K to populate the excited states, the 
ESR data for the various R impurities can be explained 
and the full set of CFPs determined. 

The last term in Eq. ([1]) is usually of secondary impor- 
tance. ESR is the second technique known to us where 
this term cannot be ignored. The other examples are the 
crystalline field potential of PrOs4Sbi2 and PrFe4Sbi 2 
measured by inelastic neutron scattering. 12, 1^, 2^- In 
those compounds the Bg-term rules out the non-Kramers 
doublet T3 as the ground state, in favor of the Fi singlet. 



EXPERIMENTAL 

Single crystals of Cei-^R:rFe4Pi2 (R = Nd, Dy, Er, 
Yb; X < 0.003) were grown in a molten Sn flux accord- 
ing to the method described in Ref. [2lj. Whithin the 
accuracy of microprobe analysis the crystals studied are 
found to be uniform. The R concentrations were deter- 
mined from the H and T-dependence of the magnetiza- 
tion, M{H,T). M{H,T) measurements were taken in a 
Quantum Design MPMS SQUID dc-magnetometer. The 
crystals used were about 2x2x2 mm^ with perfect nat- 
ural crystallographic grown faces. The cubic structure 
(space group Im3) and phase purity were checked by 
x-ray powder diffraction. The ESR spectra were taken 
in Bruker X (9.48 GHz) and Q (34.4 GHz) band spec- 
trometers using appropriated resonators coupled to a T- 
controUer of a helium gas flux system for 4.2 ^ T < 300 
K. The R^+ resonances show dysonian (metallic) line- 
shape {A/B « 2.5) corresponding to a microwave skin 
depth {5 = l/(7r/io(Ti')^/^) smaller that the size of the 
crystals. (2^ The low-T metallic character of the com- 
pound is associated to the thermally activated conduc- 
tivity (~ 10"'^ (ficm)^^) reported for this material at 
low-T. [21|] 



RESULTS 

Fig. H] shows the T-dependence of the X-Band ESR 
linewidth, AH, for the Kramers doublet ground state 
of Dy'^^ and Er'^^ in CeFe4Pi2. Within the experimen- 
tal accuracy, the linear T-term is negligible at low-T, 
in agreement with previous measurements. fl7| This in- 
dicates that there is no spin-lattice relaxation via an 
exchange interaction with the conduction-electrons (ce) 



(Korringa relaxation). [231 124| Q-Band data (not shown 



here) are similar to the data presented in Fig. [T] with 
slightly larger (< 15%) residual linewidth, AH{T = K), 
i.e., no inhomogeneous broadening is observed. Thus, the 
exponential increase of AH at high-T results from an ho- 
mogeneous line broadening due to a phonon spin-lattice 
relaxation process involving the excited CF levels (see 



below). [li, [ill For Yh^+ in CeFe4Pi2 a T-independent 
(not shown) resonance of AH = 8(2) Oe corresponding 
to a Kramers doublet ground state was observed up to 
T ~ 40 K. Fig. [2]displays M {H, T) for the same samples. 
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FIG. 1: (color online) T-dependence of AH of X-band ESR 
for Dy^+ and Er^+ in Cei_:i,RxFe4Pi2 (R = Dy, Er). The solid 
lines are fits to Eq. ([6]) leading to the following parameters 
for Dy: a = 4.0(4) Oe, Ai = 40(8) K, A2 = 135(30) K, ci 
= 0.0015(2) Oe/K^ ca = 0.0020(2) Oe/K^ and for Er: a = 
3.3(3) Oe, Ai = 85 (15) K, Aa = 300(100) K, ci = 0.0003(1) 
Oe/K^ C2 = 0.0002(1) Oe/K^ 

For Nd^+ in CeFe4Pi2 the ground state corresponds to 
an anisotropic quadruplet. The value anisotropy has 
been obtained by measuring the two allowed transitions 
within this quadruplet for the field in the (110) plane at 
T = 4.2 K [13. 

The resonances associated to the above ESR data cor- 
respond to the R-^+J = isotopes. We have also ob- 
served the resonances corresponding to various R^+ iso- 
topes with / ^ 0, that, at low-T, show the same features 
already reported. |17lll8l| Furthermore, the T-dependence 
of the ESR intensity for the observed resonances follows 
approximately a Curie- Weiss law at low-T. This indi- 
cates that the resonances arise from the ground state of 
the CF split J-multiplet. The measured 5- values and 
degeneracy of the ground states are displayed in Table I. 



ANALYSIS AND DISCUSSION 

We now add the Zeeman term gj/^sH • J to Eq. ([T]), 
where gj is the Lande 3-f actor, ^lb the Bohr magneton, 
J the total angular momentum for each R ion and H 
is the c?c-magnetic field. Following Lea, Leask and Wolf 
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TABLE I: ESR and CFPs for Cei_3:R^Fe4Pi2 (R = Nd, Dy, Er, Yb). Ground state degeneracy is abbreviated "gsd" and 
"Anis" denotes anisotropic ground state. The symbol (*) denotes a result obtained from Dy'^"'' data. 
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FIG. 



H(T) 



M{H,T) for Dy='+, Er^+ and Yb='+ in 



Cei-a;RiFe4Pi2 (R = Dy, Er, Yb). The dashed curves are 
the calculated M (H, T) from Eq. (O using the CFPs from 
Table I. The contribution of the sample holder and host lattice 
to the measured magnetization has been subtracted. 



(LLW) [25| the Hamiltonian can be parametrized as 



He 



FZ 



M^|(i-|y|) 

+gjfiB'il ■ J, 



(2) 



where we denoted O? - 



- 50| by 0% and similarly the sixth 



order terms by 0§ and Og, respectively. The coefficients 
of Eq. (dl) are rewritten as B| = (1 - |y|)xVF/F^, = 
(1 - |y|)(l - |x|)l^/Fg and Bi=yW /Yl. The coefficients 
F™ are tabulated in Ref. [26| for various values of J. 
The above is a generalization of the LLW Hamiltonian 
that includes the Og-term.[l, [2^ Our parametrization 
is slightly different from that in Ref. ^ and has the 
advantage that the entire range of the CFPs is accounted 
for within the finite intervals (-1 < x < 1) and (-1 < y < 
!)■ 



By diagonalizing He fz we obtained the CF wave func- 
tions and energies for each of the R in units of as a 
function of x and y. Then, for a small H the doublet 
ground state (F^, i = 5,6 or 7) the g- value can be cal- 
culated (g — 2gj\(yi\S:!,^i)\). For finite field and at res- 
onance g can be obtained from the Zeeman splitting of 
the doublet, /^E{H) — hv — gfiBH. Fig. [3] shows the x 
and y dependence of the g-value for the ground state of 
Er3+ (J = 15/2; g., = 6/5; > 0) in a color scale. For 
y = and variable x, we obtain the expected g-value of 
6.000 (orange) for the Fg and 6.800 (yellow) for the Fy 
doublets. [251 The white region in Fig. [3] corresponds to 
the quadruplet ground states. Thus, a g-value of 6.408 
corresponds neither to a Fg nor to a Fy (y = 0). For y ^ 
the g-value decreases and approaches to zero for large 
values of y (black region). Hence, the measured g-value 
of ^ 6.4 for Er'^+ corresponds to a doublet ground state 
with y ^ 0. Such a g-value is obtained for the set of (x,y) 
values indicated by the dashed blue line in Fig. [31 The 
results shown in Fig. [3] do not depend on the sign of y. 




FIG. 3: (color online)The color scale shows the ground state 
theoretical g-values for Er^+ (J = 15/2; gj = 6/5, W > 0) as 
a function of (x,y). The blue dashed line indicates the set of 
(x,y)-values corresponding to the experimental g = 6.40. The 
red line correspond to the (x,y)-values for Dy^"^ (J = 15/2; 
gj = A/Z, W > 0) and measured g = 7.438 (the experimental 
uncertainty of the t;- value is about the width of those lines). 

The same procedure was followed with the measured 
g-values for Dy3+ (J = 15/2; gj = 4/3) and Yh^+ (J = 
7/2; gj = 8/7) impurities. To be able to present the 
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(/-values of Dy^'^ in Fig. 3, we re-scaled them by the gj 
ratio between Dy^+ and Er3+ (gj^/gf'' = 10/9). The 
(x,y) values corresponding to the experimental 51-value 
for Dy^+ are given by the red curve in Fig. \3\ The 
results for Yb^+ are similar, but not shown in Fig. [31 

The parameters are angular momentum effective 
values of the actual CFP defined in real space. The 
and A™ are related by = (r")6'„A™.[2i] Here 6'„ 
is a geometrical factor arising from the addition of angu- 
lar momenta. The substitution of a weakly intermediate 
valence Ce ion by a R'^"'" impurity may distort the elec- 
tron density in the neighborhood of the defect. However, 
the host perturbation by the R'^"'" impurities should be 
comparable for Er'^^, Dy'^^ and Yb^"*" ions, therefore, the 
actual CFPs should not depend much on the ion R. 
Thus, if Ri and R2 denote two rare earth impurities, it 
is possible to relate their CFPs ^ 
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we obtain the following relations among the sets of pa- 
rameters (x2,y2,W2) and (xi,yi,Wi) for the two ions 
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and 
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where ^, 7 and 77 are geometrical parameters that only 
depend on On and (their values are tabulated in Ref. 
[26|). On the other hand, /3 and S depend on the expecta- 
tion values (/""(R)). (3 enters the expression for X2 and y2, 
while in order to obtain W2 also 6 is needed. The values 
for (r"(R)) have been computed in Ref. for the free 
(unperturbed) rare earth ions. In general, the (r"(R)) 
values depend on the host, in particular wether it is an 



insulating [23| or a metallic [29| environment. Their val- 
ues may be obtained from ab-initio calculations, which 
are beyond of the scope of this work. Nonetheless, /3 
and S depend on the (r"(Ri))/(r"(R2)) ratios that, for 
Dy"^"^, Er^+ and Yb'^+, in a given lattice (whether insu- 
lator 27, [2^ or metals [2^) present differences smaller 
than 5%. In other words, the changes of (r"(R)) from 
the free ion values are about the same for the various R 
when located in the same environment. As CeFe4Pi2 is 
a small gap semiconductor, for f3 and 6 we shall assume 
values close to those for an insulator. Here we assume 
that the values are within ±10% of the insulating ones. 

M 

The blue curve in Fig. d] again shows the (x,y) pa- 
rameters for Er^^^ (see Fig. [3]). Using Eq. [H the set 
of (x,y) values for Dy'^^ and Yb'^^ that satisfy the mea- 
sured ground state g-values may be transformed to the 
(x,y) space corresponding to Er^+. The results for Yb'^+ 
and Dy'^+ are shown in Fig. [4] by the red and black lines, 
respectively. The width of these lines includes the un- 
certainty of P and experimental error bars of the mea- 
sured (7- values. Notice that the lines for Dy''+, Er-^^ 
and Yb'^+ all intersect at a single point (x « —0.16(3), 
y « 0.45(3)). The three ions have the same charge 
and a similar size, therefore, we may assume that the 
actual CFPs are about the same for these impurities 
in CeFe4Pi2. This suggests that the ratios involving 
the actual CFPs are A|(r4(R))/^§(r6(R)) « -2.0 and 
Aq/Aq « 10 for the three impurities. Now these Er^^^ 
(x,y)-values are transformed back to obtain the (x,y)- 
values for Dy'^^ and Yb'^+, which are listed in Table I. 
Notice that the (x,y)-values for Dy'^+ and Yb'^"'" are ob- 
tained by using their experimental (7-values and the as- 
sumption of similar crystal fields. 
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FIG. 4: (color online) Set of (x,y)-values, satisfying the 
ground state g-values of the studied R^"*" ions, transformed 
into the Er^^ (x,y)-space by eq. 4. The open circle indicates 
the point where Dy^^, Er^^ and Yb'^"'" share the same ratios, 
A'4{r*{R))/A^e{r'^{R)) ^ -2.0 and Al/A^^ « 10 (see text). 



The above is valid independently of the energy scaling 
parameter W. The values of W for all three impurities 
can be determined if the W for one of them is known (see 
Eq. (O). W can be estimated from the T-dependence of 
AH data. The solid lines seen in Fig. [1] are the best fit 
of the measured AH for Dy^+ and Er^+ in CeFe4Pi2 to 
the expression 
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where a is the residual linewidth. The relaxation is 
through phonon modes and requires the coupling of 
phonons between the ground and excited CP states. We 
consider here the two lowest excited CF states with non- 
vanishing matrix elements from the ground state and de- 
note the excitation energies with Ai,2, respectively. The 
coefficients ci,2 are given by {3kl/2Trh^pv-')Ml2 {p is the 
host density, v the sound velocity and 2 the sum of 
the square of the matrix elements of the dynamic-crystal 
field potential) . [3, The parameter values resulting 
from the fits are given in the caption of Fig. [1] 

By using the values of (x,y), Ai and A2 (see Fig. [T]) 
we obtain Woy = 0.92(16) K and Wet = 1-6(3) K. The 
resulting energy levels for Dy^+ and Er'^+ ions are shown 
in Fig. O On the other hand, using Woy and Eq. [5] we 
can obtain W for the other ions. In particular, for Er^^^ 
we obtain W^^. — 1.3(4) K, where the error bar includes 
20% of experimental errors and 10% from the uncertainty 
of S. We see that, within the error bars, the values Wet 
and W^^ agree. Using again Woy we determined Wyh = 
7(2) K. The Yb'^+ energy levels are also shown in Fig. [5l 
Once the set of (x,y) and W parameters are known for 
a given R, their corresponding CFFs are calculated 
(see Table I). 

Assuming that the (r"(R)) values are, within 10%, of 
their values in insulators, [28| the actual CFFs A™ can 
be estimated 



AI ^ -33(10)K/a4, 
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(7) 



where ao is the Bohr radius. The main sources of error 
are the experimental Ai and theoretical S uncertainties. 

For Nd^+ the (x,y) values that account for the mea- 
sured anisotropic (/-values of the quadruplet ground state 
are given in Table L These values correspond to a point 
in the Er'^+ (x,y) space (filled circle in Fig. |4]), which is 
different than that for the other R'^^ ions. This suggests 
that the large Nd^+ ionic radius, as compared with those 
of the other R'^+ and the intermediate valence of the Ce 
ions, probably causes a large local crystal distortion close 
to the Nd^+ site. Additional experimental information, 
involving the excited CF levels, would be needed to de- 
termine WncI in this compound. Therefore, the complete 
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FIG. 5: (color online) CF energy levels for the studied R'^"'" 
ions in CeFe4Pi2. For Nd^"*" we used arbitrarily WncI = -5 
K. Thus, if the actual VFjvd value is known the Nd^^ energy 
levels should be scaled by -WNd/5 K. The high of the dotted 
line boxes indicate the uncertainty of the energy levels. 



set of actual CFFs A™ for Nd'^^ cannot be given. Nev- 
ertheless, the (x,y) values for Nd'^+ are compatible with 
WNd < 0. In Fig. [5] we present the Nd'^"'" energy levels 
using arbitrarily WNd = -5 K (see caption of Fig. (5]). 

The CF splittings of the J-multiplet determine 
M{H,T) 



M{H,T) 



^,^+V,(H)e-^-(^)/fe^ 

^2^+1 e-EdH)/kBT 



(8) 



where mi{H) and Ei{H) are the magnetization and en- 
ergy eigenvalue of each eigenstate of the Hamiltonian 
(Eq. m computed at a, finite H using Eq. [5] and the CFFs 
B™ given in Table I. The dashed curves in Fig. [2] show 
the calculated magnetization, M{H, 2 K) and M{H, 12 
K), for two concentrations of Dy'^+, Er'^+ and Yb'^+ as 
compared to the experimental data. In all cases the sam- 
ple holder diamagnetism was previously determined and 
subtracted from the total magnetization. Also, the para- 
magnetic contribution of the undoped CeFe4Pi2 host lat- 
tice was measured and subtracted. At low-T (2 K ^ T 
^ 12 K) that magnetization is less than 10 % of the sam- 
ples doped with Er {x = 0.0017) and Dy {x = 0.0034) 
and ~ 30 % of the one doped with Yb {x = 0.0023). 
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CONCLUSIONS 



The filled skutterudite CeFe4Pi2 compound is a small 



gap (~ 1500 K) semiconductor. [2l| Hence, the R + spin- 



lattice relaxation via an exchange interaction with ce is 
inhibited (Korringa process), 23, 24] since the ce must be 
promoted via exponential activation. This is verified by 
the absence of a linear T-term in our low-T AH data 
(see Fig. [1]). Similarly, a 5-shift (Knight shift) 24] is 
not expected. Therefore, the shift of the g-value of the 
Kramers ground-doublet relative to that in Oh symmetry 
(y = 0) is due to the B^(0| - 0^)-term in Hcf- For 
symmetry the Kramers doublet g-values are unique (in- 
dependent of the CFPs) [iij and the exchange coupling in 
a metallic host is simple obtained from the g-shift of the 
resonance. Our calculation showed that the presence of 
the new term results always in an isotropic 5-value and a 
negative g-shift for the doublet ground states for J — 7/2 
and 15/2. For impurities in metallic hosts with T^ sym- 
metry, when studied by ESR, a negative g-shift results in 
a complication to evaluate the sign and magnitude of the 
exchange interaction between the R'^^ localized magnetic 
moment and the ce. 

In summary, in this work we measured the ESR for 
Dy'^+, Er'^+ and Yb'^+ ions doped into the filled skut- 
terudite CeFe4Pi2 with T/i structure. We obtained the 
three CFPs B™, determined the CF ground state, ex- 
plained the unexpected Er'^+ g-value, and found the CF 
overall splitting for the J-ground state multiplet. With 
the obtained CFPs we could fit the low-T M {H, T) of 
the crystals used in the ESR experiments. Moreover, our 
working assumption that the actual CFPs A™ are about 
the same for Dy^^+j Er^+ and Yb'^+ in this compound 
turned out to be very plausible. A similar work could be 
carried out on undoped compounds such as LnFc4Pi2, for 
Ln = Nd, Gd, Dy, etc., all Kramers ions with magnetic 
ground multiplet. Besides, this work and our preliminary 
ESR data in the doped unfilled skutterudites CoSbs put 
in evidence the importance of the extra B\{0\ — Og)- 
term in Hcf for compounds with T/j symmetry. In ad- 
dition, we emphasized the extra caution we need to have 
when ESR is used to determine the exchange parameter 
in metallic compounds with T and T/j symmetry. 
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